First order phase transition from the vortex liquid to an amorphous solid 
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We present a systematic study of the topology of the vortex solid phase in superconducting 
Bi2Sr2CaCu2 08 samples with low doses of columnar defects. A new state of vortex matter im- 
posed by the presence of geometrical contours associated with the random distribution of columns 
is found. The results show that the first order liquid-solid transition in this vortex matter does not 
require a structural symmetry change. 



As a consequence of the interaction between vortices 
and quenched disorder in superconducting samples, the 
vortex liquid solidifies through a first order phase transi- 
tion (FOT) [1-3] into a topologically ordered solid with 
elastic deformations, the Bragg glass phase [4,5]. Arti- 
ficially introduced correlated defects, as those produced 
by heavy ion irradiation, transform the first order liquid- 
solid transition into a second order one in which the solid 
vortex phase, the Bose Glass [6], has no topological quasi 
long range order. However, studies [7] of the phase dia- 
gram of Bi 2 Sr 2 CaCu20 8 (BSCCO-2212) have shown that 
the FOT is robust in the presence of low density of colum- 
nar defects. In this case the FOT is preserved even for 
vortex densities comparable to those of defects. 

High dose of columnar defects has two main effects on 
the vortex system. On one hand, its random distribution 
suppresses the FOT; on the other hand, the correlated 
potential hardens the vortex lines inducing a divergent 
vortex tilting elastic constant, C44, at the Bose Glass 
transition. The results in Ref. [7] indicate that the dis- 
order induced by a low density of columnar defects is 
too weak to change the order of the liquid-solid transi- 
tion. This could imply that the quasi long range order 
of the Bragg Glass is preserved despite the presence of 
correlated defects. 

The aim of this work is to establish the correlation be- 
tween the space configuration of columnar defects and 
the vortex structure created in the presence of those de- 
fects. This goal is achieved visualizing the vortex struc- 
ture by means of the magnetic decoration technique in 
samples where the liquid-solid phase transition has been 
characterized by differential magneto optics (DMO) [8]. 
The observed vortex topology is discussed in terms of the 
nature of the liquid-solid transition. 

The BSCCO-2212 samples used in this work were ir- 
radiated at GANIL by 1 GeV Pb ions parallel to the 
c axis of the crystal. In order to compare the vortex 
structure in the presence of columnar defects with that 



of the pristine sample, the irradiation was made through 
an ad-hoc stainless steel mask that allows to have irradi- 
ated and non irradiated regions in the same sample (see 
inset of Fig 1(a)). We have investigated four samples ir- 
radiated with doses corresponding to B$ = 5, 10, 50 and 
100 G, where B&, the matching field, is usually defined 
as n co ;$o with n co i the density of columns and $0 the 
flux quantum. The H — T phase diagram obtained by 
DMO imaging [9] for these samples is shown in Fig. 1(a). 
It is important to remark that this high sensitivity tech- 
nique shows that the FOT is not only preserved in the 
presence of low density of columnar defects but the tran- 
sition temperature is raised when compared to that of 
pristine samples. For samples with B$, — 5 and 10 G the 
melting line is first order and for the J3$ = 100 G sample 
is continuous in the range of fields investigated by mag- 
netic decoration. For the B$> = 50 G sample a critical 
point at 70 G separating the FOT at lower fields from a 
continuous one at higher fields is found. The magnetic 
decoration was performed after cooling the sample in the 
presence of a field, FC experiment, down to 4.1 K. The 
results obtained in the non irradiated regions confirm the 
high quality of the crystals. 

Figure 1 (b) shows the vortex structure in the non irra- 
diated and irradiated regions obtained by magnetic dec- 
oration after FC the Bq, = 10 G sample in a field of 40 G. 
A sharp boundary separates the Bragg Glass in the non 
irradiated region from a structure with no long range 
crystalline order in the irradiated one. 

The experimental results show that the vortex struc- 
ture is polycrystalline for B > £>$ and amorphous in the 
range B < _B$ . Examples of Dclaunay triangulations of 
polycrystalline vortex structures observed by magnetic 
decoration after FC the Bg, = 10 G sample in fields of 
80 G and 40 G are shown in Fig. 2(a) and (b), respec- 
tively. Fig. 2(c) shows the vortex structure for B = £>$. 
The black dots indicate topological defects (non-sixfold 
coordinated vortices) . Figure 3 depicts decoration images 
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of amorphous structures obtained in the regime B < _B$; 

(a) after FC the i?$ = 50 G sample in a field of 30 G and 

(b) after FC the i?$ = 100 G sample in a field of 80 G. 
Contrary to what is often observed in homogeneous 

materials, we argue that the nature of the polycrystallinc 
structure shown in Fig. 2 is an intrinsic property of the 
vortex matter in the presence of randomly distributed 
correlated defects. To demonstrate this we have used 
the dynamic annealing method [10,11] in which the mag- 
netic field is tilted back and forth before the decoration 
is made. This procedure is known to remove the grain 
boundaries in NbSe2 and to reorder the crystalline vortex 
lattice [10] in non irradiated BSCCO-2212 samples. In 
contrast, we find that the polycrystalline structure of the 
irradiated region is not altered after applying this anneal- 
ing procedure. This supports that the observed structure 
is the equilibrium one rather than a mctastable configu- 
ration resulting from the typical nuclcation and growth 
processes. 

The above discussion suggests that the grain bound- 
aries of the polycrystalline structure are fixed in space 
by the random distribution of columnar defects. If this 
were the case the area of the grains for a given £?$ would 
be field independent contrary to what is observed in ho- 
mogeneous low T c superconductors. To verify this, we 
studied the statistical distribution of the grain areas as 
a function of the applied field for B$ = 5 and 10 G sam- 
ples. The grain boundaries are defined as the interface 
between crystallites which orientations differ by 10 de- 
grees or more. The areas of the different grains are mea- 
sured counting the number of vortices within them. The 
analysis of many pictures shows that the histogram of the 
grain sizes is broad with upper and lower limiting values 
for the grain areas for each applied field. The number 
of vortices, N vg , within the largest and smallest grains 
(upper and lower limits of the histogram) are propor- 
tional to B as shown in Fig. 4 for B& = 10 G sample. 
Identical behavior is found in the £?$ = 5 G sample with 
the shift of the distribution of grain areas corresponding 
to the lower irradiation dose. Furthermore, the average 
number of vortices per grain (defined as the ratio be- 
tween the total number of vortices within grains and the 
number of grains) is also found to be proportional to B, 
see Fig. 4 (solid squares). These results make evident 
that the size of the different grains for a given B$ is field 
independent contrary to what is found in low T c super- 
conductors, where the average grain size grows with the 
magnetic field [12]. The invariance of the grain size with 
B in the irradiated regions strongly supports that the 
grain size as well as the space distribution of the grains 
are determined by the landscape of the columnar defects. 

In order to understand the relation between grain 
boundaries and the distribution of correlated defects we 
have generated a set of random points simulating the po- 
sitions of columnar defects, see Fig. 2(d) (white dots). 
This configuration is found to be equivalent to that of 



columnar defects detected after etching mica irradiated 
by heavy ions with comparable doses, see inset in Fig. 
4. It is interesting to remark that the distribution of 
columnar defects is not homogeneous on the scale of the 
lattice parameter in the range of fields investigated. The 
in-plane inhomogeneous distribution and the correlated 
character of the columnar defects suggest the presence of 
a network of contours extending throughout the thick- 
ness of the sample associated with the localization of 
vortices on columns. In Fig. 2(d) the contours are repre- 
sented by the gray connected areas indicating the range 
of interaction ( penetration depth) of vortices localized 
on defects. These vortices create energy barriers that 
inhibit the propagation of the crystalline vortex lattice, 
giving rise to the polycrystalline structure observed in the 
regime B > _B$. Thus, the observed grain boundaries in 
samples with low irradiation dose have a different origin 
than those associated with metaestable states in homoge- 
neous systems [11]. The apparent grain boundaries in the 
irradiated samples are field independent vortex contours 
fixed in space by the random distribution of columns. 
This prompts us to propose a new state of vortex mat- 
ter where the total number of vortices is divided in two 
species, one corresponds to the fraction of vortices local- 
ized on contours, p CO nt, and the other forming the crys- 
tallites surrounded by the contours, p crys - 

Based on the above discussion we present a simpli- 
fied picture that allows to predict the dependence of 
the number of topological defects (non sixfold coordi- 
nated vortices) of the polycrystalline structure, Na e f, on 
B and The total number of defects is counted in 

an arbitrary large area that contains many grains. The 
reported data was obtained from an inspected area of 
the order of 10 4 \x m 2 , much larger than the maximum 
grain size for the different B$ investigated. From the 
inspection of many images of the vortex structure in 
the irradiated regions we find that most of the topo- 
logical defects are located along the contours separat- 
ing the crystallites. Thus, as a first order approxima- 
tion we assume that Ndef is proportional to the num- 
ber of vortices localized on the contours, N cont . The 
total length of the contours within the inspected area 
can be written as L — N cont cto where ao is the aver- 
age vortex distance. The experimental results showing 
the field independence of the areas of the crystallites (see 
Fig. 4) imply that L is only function of B$. There- 
fore N de f oc N cont oc L(B^)/a oc L(B^)B X I 2 . On the 
other hand, the self similarity of the space distribution 
of random columnar defects imposes that the areas en- 
closed by the contours scale with and consequently 
L(B$) oc B§ 2 . Therefore, the fraction of vortices in- 
volved in topological defects, p de f, is 

Pdef = N def /N v cx {B^/B) 1 ' 2 (1) 

In Fig. 5 we plot pdef obtained from Delaunay triangu- 
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lations of vortex structures for all B and investigated. 
The black line corresponds to the square root dependence 
predicted by Eq. 1. Despite the crude approximations, 

1/2 

Ndef oc N cont and L(£>$) cx , we see that for < B 
the experimental results are well described by this simple 
model. On the other hand, the deviation of the exper- 
imental data from the dependence predicted for pdef in 
the range B^jB > 1 is evident. 

The failure of the model to describe the amorphous 
state is expected, since for fields B < B$> a significative 
fraction of the areas enclosed by the contours fixed in 
space are not large enough to allocate vortices with crys- 
talline symmetry, p crys ~ 0. In this limit, despite most 
of the vortices remain pinned on columns, p con t ~ 1, 
the number of topological defects is not proportional to 
L(£>$)/a ; i.e. a is not an appropriate length scale to 
apply the scaling argument. 

The B<t>/B independence of the fraction of defects in 
the amorphous state (pdef ~ 0.5) is remarkable, see Fig. 
5. To investigate the origin of this result we have gener- 
ated a random distribution of points with the constrain 
that the distance between near neighbors is not less than 
50% of the lattice parameter of a perfect lattice with the 
same density. The constrain is chosen to allow the fluctu- 
ations of the average distance between points to be that 
observed in the amorphous vortex structure. It is found 
that the fraction of defects obtained from Delaunay tri- 
angulation of this simulated distribution of points is the 
same as that for the vortex structure for B < S$, indi- 
cating that in this regime most of the vortices are pinned 
on a random distribution of columnar defects. 

The analysis of the vortex structure for B > i?$ as 
well as for B < £>$, Fig. 5, strongly supports that in 
both limits the topological defects are associated with 
vortices localized on columns. For B > the topo- 
logical defects of the vortex structure are mainly due to 
plastic distortions induced by the presence of crystallites. 
In the other limit, B < i3$, the space within the contours 
does not allow to form crystallites, p cry s ~ 0, and the de- 
fects are associated with vortices localized on a fraction of 
columns compatible with the constrain mentioned above. 
It is clear that B/B& « 1 marks the transition from a 
topological state described by a collection of crystallites 
to a state characterized by an amorphous structure, both 
determined by the random distribution of columns. 

As a consequence of the results discussed above we 
see that in the limit B < _B$ no crystalline structure 
can be accommodated within the small contours deter- 
mined by the columnar defects, p cr ys ~ 0. Thus, an 
amorphous vortex structure characterizes the solid state. 
This is particularly relevant when analyzing the interre- 
lation between the solid structure and the nature of the 
vortex liquid-solid transition. The critical point at 70 G 
in the phase diagram of the B$ = 50 G sample indicates 
that most of the field range where the first order melting 



takes place corresponds to the regime B < Thus, 
the observed amorphous vortex structure, see Fig. 3(a), 
rules out the widely accepted correlation between a first 
order liquid-solid transition and a structural symmetry 
transformation. Moreover, we observe no difference be- 
tween the topological order in Figs. 3 (a) and (b) despite 
the fact that the former structure has solidified through 
a FOT whereas the latter through a continuous phase 
transition. 

In summary, a new vortex matter with two types of 
vortices has been discovered, giving support to the re- 
cently suggested vortex porous structure [9] and the in- 
terstitial liquid [13]. One type is associated with vortices 
localized upon the topological contours formed by the 
random distribution of columnar defects and the other 
with vortices forming the crystallites. The relative distri- 
bution of vortices in these two classes explains the topo- 
logical transformation of the solid from a Bragg Glass 
= 0) into a polycrystal < B) and then into 
an amorphous structure (B$ > B) as the density of 
columnar defects is increased. The lack of correlation 
between the structure of the vortex solid and the order 
of the melting transition opens an important question 
on the microscopic mechanism that triggers melting in 
vortex solids. The detection of a FOT between an amor- 
phous solid and a liquid makes evident that the first order 
melting transition in vortex matter does not require the 
widely accepted long range structural order of the solid 
state. 
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FIG. 1. (a) Phase diagram obtained by DMO for the pris- 
tine BSCCO-2212 sample and for samples with columnar de- 
fect densities corresponding to B<s> = 5, 10, 50 and 100 G. Solid 
lines indicate first order liquid-solid transition while dotted 
lines correspond to continuous ones. Inset: Schematics of 
the sample with the mask used for irradiation, (b) Mag- 
netic decoration at 40 G in the sample with = 10 G. The 
dashed line shows the boundary between the irradiated and 
the non-irradiated regions. 

FIG. 2. Delaunay triangulations of the vortex structure in 
a sample with B& = 10 G at fields of: (a) 80 G, (b) 40G 
and (c) 10 G. The black dots indicate the non sixfold coordi- 
nated vortices. The gray lines in (a) and (b) depict the grain 
boundaries, (d) Set of random points (white dots) generated 
numerically simulating a distribution of columnar defects with 
B$ = 10 G. The gray areas indicate the range of interaction 
of vortices localized on defects, see text. 



FIG. 3. (a) Structure of the vortex solid after a FOT at 
B = 30 G in the B<s> — 50 G sample, (b) Structure of the 
vortex solid after a continuous liquid-solid transition at 80 G 
in the B$ = 100 G sample. 

FIG. 4. Number of vortices within the largest (triangles), 
the average (squares) and the smallest (circles) grains as a 
function of B for B<s> — 10 G. Inset: Columnar defect distri- 
bution in mica irradiated with heavy ions. 

FIG. 5. Fraction of defects, pdej as a function of B&/B 
for all Bi> investigated. The black line corresponds to 
pde} oc {B^/B) 1 ^ 2 and the gray dotted line to p de f — 0.5. 
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Figure 1, M. Menghini etal. 




Figure 3, M. Menghini etal. 
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Figure 4, M. Menghini etal. 




Figure 5, M. Menghini etal. 



